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ABSTRACT 

We present mid-infrared spectra of sixteen optically faint sources with 70/im 
fluxes in the range 19mJy < f„ (70/im) < 38m Jy. The sample spans a redshift 
range of 0.35 < z < 1.9, with most lying between 0.8 < z < 1.6, and has in- 
frared luminosities of 10 12 — 1O 13 L . Ten of 16 objects show prominent polycyclic 
aromatic hydrocarbon (PAH) emission features; four of 16 show weak PAHs and 
strong silicate absorption, and two objects have no discernable spectral features. 
Compared to samples with f„(24//m) > lOmJy, the 70/im sample has steeper IR 
continua and higher luminosities. The PAH dominated sources are among the 
brightest starbursts seen at any redshift, and reside in a redshift range where 
other selection methods turn up relatively few sources. The absorbed sources 
are at higher redshifts and have higher luminosities than the PAH dominated 
sources, and may show weaker luminosity evolution. We conclude that a 70/im 
selection extending to ~ 20mJy, in combination with selections at mid-IR and 
far-IR wavelengths, is necessary to obtain a complete picture of the evolution of 
IR-luminous galaxies over < z < 2. 
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Introduction 



Among the most important cosmological results of the last few decades was the dis- 



covery by th e Cosmic Background Explorer 



wavelengths (IPuget et al.lll996l ; lHauser et al. 



COBE) of a background radiation at infrared 



19981 ). This background is comparable in in- 



tensity to the integrated optical light from the galaxies in the Hubble Deep Field, implying 
that the star formation rate density at z > 1 was more than an order of mag nitude higher 



than l ocally, and that most of this star formation wa s obscured. Later surveys (lAussel et al. 



1999 



Dole et al.l l200ll ; iRowan-Robinson et ajj l2004j ) resolved the bulk of this background 
into a population of distant IR- luminous galaxies (LIR Gs, L ir > 10 n L M ) which undergo 



stron g luminosity evolution with redshift ((1 + z)~ 4 , e.g. iPozzi et al.l 12004 ; iLe Floc'h et al. 



20051 ). re aching a comovin g density at least 40 times greater at z ~ 1 than in the local 



Unive rse (lElbaz et al.ll2002f ) . Rey iews of their properties can be found in lSanders Mirabel 
Jl996h and lLonsdale et al.l J2006h . 



Significant effort has been devoted to understanding the mechanis ms driving the evo- 



lution of LIRGs. At low redshift, they are almost invaria bly mergers (ISurace et al.l 11998 



Farrah et al.l 



2001: 



Bus house et alll2002l; IVeilleux et al.ll2002T). powered mainly by star f orma- 



tion ( Genzel et al 19981 ; Rigopoulou et al. 1999 ; Imanishi et al. 2007 ; Vega et al. 2008 ). and 



reside in average density environments ( Zauderer et al. 2007h . LIRGs at high redshift also ap 



pear to be mainly starburst dominated merging systems (IFarrah et al. 



2003b; Smail et al 



2007 



Bridge et al. 



20041; iTakata et all2006l : lBorvs et alfeood : 



2002; 



Valiante et al 



Chapman et al. 



2007 



Bertaetal 



20071 ) , though there are signs of differences compared to their low red- 



shift counterparts; for exa mple, weak X-ray emission (jFranceschini et al. 



20031 ; llwasawa et al.ll2005l ). different modes of star formation (IFarrah et al 



2003 



dency to reside in overdense regions (IBlain et alll2004l ; IFarrah et al 
20071 ). 



2006 



Wilman et al. 



2008 ), and aten- 



Magliocchetti et al. 



Controversies remain, therefore, over how LIRGs may or may not evolve with redshift. 
Part of the reason for this is that an efficient census of LIRGs at z > 0.5 is difficult, as surveys 
conducted in a single IR band can miss a significant fraction of the LIRG po pulation. For 



exam 


pie, sub-mm surveys find large numbers of obscured starbursts at z > 1 ( 


Barger et al. 


1999; 


Chapman et al. 


2005; 


Aretxaga et al. 


2007; 


Clements et al. 


2008; 


Dye et ; 


il. 


2008), but 



few sources at z < 1, and virtually no sources with 'hot' dust (IBlain et al.l 120021 ) . It is 
therefore essential that we survey for LIRGs in every IR band available to us and, having 
found them, systematically study them further. 



The Spitzer space telescope (jWerner et all |2004 ISoifer et al.l 120081 ) has the capacity to 
revol utionize our understanding of LIRGs. The Infrared Array Ca mera (IRAC. iFazio et al 
20041 ) and the Multiband Imaging Photometer for Spitzer (MIPS, iRieke et all 120041 ) imag- 
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ing instruments, and the Infrared Spectrograph (IRS. iHouck et al.ll2004j ) all offer dramatic 
improvements in sensitivity and resolution over previously available facilities. In particular, 
the MIPS 70/im channel is ideal for studying high redshift LIRGs; at z < 2 the rest-frame 
emission is always longward of 18/xm, giving good sensitivity to both starbursts and AGN. 
The optically faint, high redshift 70/im Spitzer sources may be a prime example of the 
sources that previous surveys in the sub-millimeter or in the infrared (IR) have missed - in 
the sub-mm because the sources harbour dust that is too hot for current generation sub- 
mm cameras to see, and in the surveys of the Infrared Space Observatory (ISO) because of 
sensitivity limits at > 20/iuQ. 



In this paper, we use the IRS to observe a sample of 16 sources selected at 70 /jm using 



data 



2003; 


Oliver et all 


2004 




Davoodi et al. 


2006alJbl: 


Shupe et al. 


2008|; 


Siana et al. 


2008|) survey. O 



2007 



luminosities and power sources seen in the optically faint 70/xm population. We assume 
Q = 1, A = 0.7 and Hq = 70 km s _1 Mpc -1 . Luminosities are quoted in units of ergs s _1 or 
of bolometric solar luminosities, where L & = 3.826 x 10 33 ergs s" 1 . 



2. Methods 

2.1. Sample Selection 

The sources are selected from the SWIRE Lockman Hole field, which covers 10.6 square 
degrees and reaches 5cr depths of 4.2/iJy at 3.6/im, 7.5/iJy at 4.5/im, 46/iJy at 5.8/im, 47/iJy 
at 8.0/im, 209/iJy at 24/im, 18mJy at 70/im, and 108mJy at 160/iin. The primary selection 
criterion for our sample is a confident detection at 70/im, so we first rejected all sources 
fainter than 19mJy at 70/im. To ensure that we could obtain mid-infrared IRS spectra 
with reasonable signal-to-noise, we also constrained the sources to have f I/ (24/im) > 0.9 mJy, 
although > 90% of sources with f 70 > 19mJy also satisfy /24 > 0.9mJy. This resulted in 
a parent sample of 1250 sources. From this, we selected optically faint sources by taking 
all sources (12 in total) with r-band magnitudes fainter than m r = 23, and including an 
additional four sources with r-band magnitudes in the range 20 < m r < 23, for a total of 16 
objects. 



^e.g. the ELAIS ISO survey reached a limiting depth of ~ lOOmJy at 90^m (jRowan-Robinson et al 

20041 ). compared to ~ 20mJy for large area surveys with MIPS at 70/im) 
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2.2. Observations 

All 16 objects were observed as part of Spitzer program 30364 with the first order of the 
short-low module (SL1; 7.4/im - 14.5/im, slit size 3.7" x 57" with 1.8" pix" 1 , R~ 60 - 127), 
and the second order of the long-low module (LL2; 14.0/xm - 21.3/xm, slit size 10.5" x 168" 
with 5.1" pix -1 , R~ 57 — 126). Eight of these objects were additionally observed with long- 
low order 1 (LL1, 19.5/im - 38.0/xm). The targets were placed in the center of each slit using 
the blue peak-up array. Each target was observed with an individual ramp time of 60s in 
SL, and 120s in LL, with the number of ramps determined by the targets 24/zm flux density. 
Details are given in Table [U 

The data were processed through the Spitzer Science Center's pipeline software (version 
15.3), which performs standard tasks such as ramp fitting and dark current subtraction, and 
produces Basic Calibrated Data (BCD) frames. Starting with these frames, we removed 
rogue pixels using the zrsc/earj^l tool and campaign-based pixel masks. The individual frames 
at each nod position were t hen combined into a single image using the SMART software 



package (IHigdon et al.ll2004h . Sky background was removed from each image by subtracting 
the image for the same object taken with the other nod position (i.e. 'nod-nod' sky sub- 
traction). One-dimensional spectra were then extracted from the images using the SPICE 
software package using 'optimal' extraction and default parameters. This procedure results 
in separate spectra for each nod and for each order. The spectra for each nod were in- 
spected; features present in only one nod were treated as artifacts and removed. The two 
nod positions were then combined. The first and last 4 pixels on the edge of each order, 
corresponding to regions of decreased sensitivity on the array, were then removed, and the 
spectra in different orders merged, to give the final spectrum for each object. 



3. Results 

The spectra are presented in Figures [1] and [2j Redshifts and fluxes are given in Table 
El and spectral measurements are given in Table [3j 



2 This tool is available from the SSC website: 



http://ssc.spitzer.caltech.edu 
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3.1. Redshifts 

We derive spectroscopic redshifts from broad emission features at 6.2/mi, 7.7/im, 8.6/xm, 
11.2/im and 12.7/im, attributed to bending and stretching modes in neutral and ionized 
Polycyclic Aromatic Hydrocarbon (PAH) molecules (the 12.7/im feature also contains a 
contribution from the [NeII]A12.81 fine-structure line), and/or a broad absorption feature 
centered at 9.7/zm arising from large silicate dust grains. For eleven sources (1, 3-7, 9, 13-16), 
an unambiguous redshift can be determined from the PAH features; the uncertainty on these 
redshifts is governed by the variations in PAH peak wavelengths seen in local galaxies and 
is of order Az = 0.02. In three cases (2, 10, 12), the PAHs are weak, and the redshifts are 
derived from a prominent silicate absorption feature. In these cases, the redshifts have a 
larger error, Az ~ 0.2, but should still be reliable. Finally, in two cases (8 & 11), no spectral 
features can be unambiguously identified; we derive tentative redshifts based on what are 
plausibly PAH or silicate features. In these cases, the error on the redshift is large, Az ~ 0.4, 
and the redshifts should be treated with caution. 



3.2. Spectral properties 



We measure the properties of the PAH features using two methods. For the 6.2/xm 
and 11.2/im PAH features, we compute fluxes and equivalent widths (EWs) by integrating 



sec 



Brandl et al. 


2006 and 


Spoon et al. 


2007) 



the continua of our sample are only weakly detected. Using this method means, however, 



(Weedman et al. 


2005; 


Brandl et al. 


2006; 


Armus et al. 


2007; 


Spoon et al. 


2007) 



7.7/im PAH feature, we do not attempt to measure EWs because of the uncertainties in 
determining a continuum baseline underneath this broad and comp lex feature. Instead , we 



measure only the flux density at the peak of the 7.7/xm feature, as in lHouck et al.l (120071 ) and 



Weedman and Houckl (120081 ) . Due to the low S/N and restricted wavelength range of the 
spectra, we cannot correct for water ice and/or aliphatic hydrocarbon absorption, although 
the effect of this lack of correction is likely to be insignificant. 

For the four objects with clear detections of the 6.2/xm and 11.2 /xm PAH features, we 
derive star formation rates via the formula from iFarrah et al.l (120071 ); this yields SFRs of 
between 100 and 30 per year . We derive star formation rates in Table 3 for all sources 
using the formula in iHouck et al.l (120071 ) : 
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logiSFfylMoyr- 1 } = log(vL v (7 '.7 >m)) -42.57 



For the four objects for which both formulae can be used, we obtain consistent results within 
the uncertainties, which are of order 50%. 

We measured the strengths of the silicate features, S s u, via: 



S s u = In ( ^ bs ^r\ ) 
\F cont {9.7nm) J 



(2) 



where F obs is the observed flux density at rest-frame 9.7/iin, and F cont is the underlying 
continuum flux density at rest-frame 9.7/im deduced fr om a spline f i t to t h e continuum on 
either side. A description of this method can be found in lSpoon et al. (b()07t ). iLevenson et al.l 
J2007h and lSirockv et all fboosh . 



Other than the PAH and silicate features, we see weak but clear detections of the H 2 S(3) 
line at 9.66/zm in two objects (6 and 14), but no other spectral features. 



3.3. SED Fitting 



We measure the IR luminosities of the sample by fitting the IR photometry simul- 
taneously with th e library of model spe ctral energy dist ributions (SEDs) for th e emission 
from a starbur st (Efstathiou et al 120001 ) and an AGN (iRowan- Robinson! Il995l ). following 
the methods in iFarrah et all (120031 ) . The fits were good, with y 2 < 2 in all cases, and the 
o bserved-frame 70/im flux gave good constraints on the IR luminosities (see also discussion 



in 



Rowan-Robinson et al. ( 20051 )). The photometry is too limited, however, to provid 



e mean- 



ingful constraints on the starburst and AGN fractions, so we only present the derived total 
IR luminosities, and not the SED fits. The luminosities are presented in Table EJ All ob- 
jects have IR luminosities exceeding ~10 11 ' 5 L Q , with most lying in the range 10 12 — 10 
making them ULIRGs. 
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4. Discussion 

4.1. Redshifts and luminosities 

The redshift distribution for our sample is shown in the left panel of Figure [3j All 
objects lie in the range 0.35 < z < 1.8. There is a peak at z ~ 0.9, a long tail up to z ~ 2, 
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and a shorter tail down to z ~ 0.4. 



The redshift range 0.5 < z < 1.5 has been a difficult one in which to select ULIRGs 
because their distances make them faint at observed-frame 10-100/zm, and the /c-correction 
that makes ULIRGs bright at observed-frame 200-1000/im does not become strong enough 
for current sub- mm cameras until z > 1.5. Moreover, this is the redshift range where the 
evolution in the ULIRG luminosity function is thought to be strongest. Therefore, our simple 
selection method, consisting of little more than a minimum 70/im flux and an optical cut, 
should prove invaluable in studying the cosmological evolution of ULIRGs from future wide- 
area surveys. It seems likely that it is the optical cut that is resulting in our sources mostly 
lying in the 0.5 < z < 1.5 range, as other spectros copic surveys of pur ely 70/im selected 
sources find a lower median redshift. For example, iHuynh et al.l (120071 ) present redshifts 
for 143 sources selected solely on the basis of 70/im flux but to a much fainter limit of 
/Vo > 2mJy, and find a median redshift of 0.64, with 79% of the sources lying at z < 1 and 
about half the sources at z < 0.5. 

Given the difficulties in obtaining spectroscopic redshifts for distant ULIRGs, it is useful 
to assess the accuracy of photometric redshifts for this type of source. In the right panel 
of Figure [31 we compar e the s pectroscopic redshifts to the photometric redshifts derived by 
Rowan-Robinson et al.l (120081 ) for those eight objects wher e the photometric redshift co de 
produced a formally acceptable fit (x 2 < 10, see discussion in iRowan- Robinson et al.l (120081 )). 
The photometric redshifts are reasonably good, given the faintness and high redshifts of the 
sample. Seven of 8 objects lie within or close to the 'catastrophic failure' boundary of 
log(l + Zphot) = log(l + z spec ) ± 0.06. Only in one case is there a clear mismatch between 
Zphot and z spec , and the photometric redshift for this object is unreliable, as it is based on 
limited data. 



4.2. Comparisons to other samples 

4.2.1. Local ULIRGs 

We first compare our mid-IR spectra to those of local ULIRGs. Our sample is faint 
at 24/xm, so detailed spectral diagnosti cs are not possible . We therefore employ a simple 
comparison using the 'Fork' diagram of ISpoon et al.l (120071 ). shown in Figure HI Our 70/im 
sam ple has similar PAH and silicate absorption properties to the 1B/1C/2B/2C classes 
from lSpoon et al.l (120071 ). This identifies our sample with moderately obscured star- forming 
sources, but not with the most heavily absorbed sources or those that contain an unabsorbed 
or silicate-emitting AGN without PAH emission. 
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4-2.2. 70fim selected samples 



The 70/im populatio n has been studied relatively little with the IRS. The only other 
published study is that of iBrand et al.l (120081 ). who select 11 sources with /70 > 32mJy and 
an r band magnitude fainter than m r = 20. Thus, the two samples make for interesting 
comparisons; our sample is 1-2 magnitudes fainter in the optical and ~ 1.6 times fainter at 
70 /im. 

The IR luminosities of both samples are comparable, with both lying mostly in the 
10 12 — 10 13 L Q range. The frac tion of sourc e s with prominent PAHs is also similar with 7 of 
11 PAH dominated sources in brand et ali ( 2008 ) and 10 of 16 in ours. There are, however, 
two areas where there are differences, albeit with the caveat of small sample sizes. The first 
is redshift distribution. The redshift distribution for the Brand et al sample is overplotted 
in the left panel of Figure [31 The Brand et al sample has a less pronounced peak, a broader 
distribution with more sources over 0.5 < z < 1 and fewer sources at z > 1.2, than does 
our sample. The second difference is the distribution of spectral types with redshift. The 
Brand et al sample shows no discernible separation of spectral type with redshift whereas 
our sample shows that all of the sources with strong PAHs, irrespective of the presence 
of silicate absorption, lie towards the lower end of the redshift range, while the strongly 
absorbed sources with negligible PAHS are all at the upper enq§- 

Both differences probably arise due to our sample reaching fainter 70 /im fluxes than the 
sample of Brand et al. In principle, surveys to fainter 70/im fluxes should include higher 
redshift, more luminous sources (see §4.2.31) . Moreover, a 70/im selection should result in 
sensitivity to different effective dust temperatures at different redshifts; at z ~ 0.7, 70 /im 
observations sample rest-frame 40/xm, while at z ~ 1.5, they sample rest-frame 28/im, so 
at z ~ 1.5 we are sensitive to sources with ~ 30K hotter dust than at z ~ 0.7. Therefore, 
higher redshift sources in a 70/im selected sample are more likely to be absorbed, AGN-like 
sources with weak PAH features, which is what we see in our sample. The absence of this 
trend in the Brand et al sample suggests that optically faint sources with f 70 > 30mJy are 
mainly ULIRGs at moderate redshift, but that at 70/im fluxes below 30m Jy we start to see 
significant numbers of heavily absorbed sources, with large masses of hot dust, at z > 1. 
Interestingly, a similar situation is seen at longer wavelengths; sub-mm surveys are adept 
at picking up sources with large masses of cold dust, but radio surveys in the same fields 
have shown that there exist pop ulations of 'hot' dust so urces at comparable redshi fts but 
with different IR spectral shapes (j Chapman et al.l (j2003al ). see also lKhan et al.l (120051 )). Our 



3 This is not a result of the 70/zm selection shifting specific spectral features in and out of the bandpass, 
as there are no prominent features that would lie at observed-frame ~ 70^m at the redshifts of our sample 
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higher redshift sources could be the lower-z tail of this radio-selected, 'hot dust' population. 



4-2.3. 24^m selected samples 



Most previous samples selected from Spitzer surveys for IRS followup use a mid-infrared 
selection based on 24/im flux. It is important, therefore, to understand whether our 70/iin 
sample differs from samples sele c ted a t 24/im. To make these comparisons, we combine our 
sample with that of iBrand et al.l (120081 ) for a total of 27 70/im selected sources, as the sample 
selections are complementary; our sample reaches fainter 70/im fluxes, but the optical cut is 
similar. 



Bright samples 

We first compare the mid-IR continuum properties of the combined 7 0/im sample to sources 



with h igh 24/im fluxes via the flux limited, f24 > 10 mJy sample in I Weedman and Houck 



( 120091 ). To perform this comparison we use the rest-frame 15/xm continuum luminosity and 
the /24//15 continuum slope. If both these rest-frame wavelengths are seen in the IRS spectra 
then we measured these quantities directly; otherwise, we estimated fluxes at one or both 
wavelengths via interpolation from a power law with a slope determined from comparison of 
observed fj, (24/im) and f„(70/im). These interpolations should be regarded with caution, as 
they are sensitive to PAH and silicate contamination of the (observed-frame) 24/im band, 
which are difficult to compute for our sample as we either lack Long-Low data, or it is of 
relatively low signal-to-noise. The ~20% uncertainty assigned in Table 3 to these interpolated 
values for the continuum slope reflects the possibility that the observed 24/im flux density 
may not be purely a measure of dust continuum emission. 

The comparison is shown in Figure [5j The continuum luminosities for the 70/im sample 
are much greater than for the 24/im sample. The median log[^L l/ (15/im)] (ergs s _1 ) for 
the 70/im sample is 44.8 compared to 43.3 for the 24/im sample. This is straightforward 
to understand. The fainter optical and 24/xm fluxes used for the 70/xm selection allow the 
discovery of IR-luminous sources to much higher redshifts so we may reasonably expect to 
see more luminous sources. 

Interestingly, however, the luminosity differences between the samples may not be as 
large for the sources with weak PAH features. For these sources, the median log(z/L^(15/im)) 
(ergs s _1 ) for the 70/xm sample is 45.4 compared to 45.0 for the 24/xm sample, and the most 
luminous sources in both samples are similar, log(i/L„(15^m)) ~ 46.2. This comparison is 
not robust, given that the weak PAH sources in the 70/im sample only have the silicate 
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feature in absorption, whereas the 24/im sample contains sources with the silicate feature in 
both absorption and emission. Nevertheless, it seems that the fainter optical and 24/im fluxes 
used for the 70/im selection do not result in discovering more luminous absorbed sources, even 
though the sources are systematically at higher redshift. We therefore suggest, with some 
reserve, that sources with weak PAHs and strong silicate absorption show weaker luminosity 
evolution with redshift than do PAH dominated sources. 

Considering the continuum slopes, we see results that would be expected from the 70/im 
selection; selecting sources at the longer wavelength favors sources with steeper spectra. For 
PAH sources, the median rest-frame ratio /24//15 for the 70/im sample is 4.5 compared to 
3.5 for the 24/im sample. For absorbed sources, the median rest-frame ratio /24//15 for the 
70/im sample is 2.7, compared to 1.7 for the 12 sources in the 24/im sample at sufficiently low 
redshifts to have a measurement. For both PAH dominated sources and absorbed sources, 
the most extreme ratios are within the 70/im sample. 

These results demonstrate a systematic difference in effective dust temperatures for the 
70/tm sample compared to the 24/im sample. For the 70/im sample, the steeper spectrum at 
rest frame ~ 24/im implies a larger dust fraction at intermediate temperatures of ~ 100 K. It 
is also notable that the PAH dominated spectra are consistently steeper than the absorbed 
spectra in both samples. In the 70/im sample, the ratio /24//15 is 4.5 and 2.7, respectively, 
and in the 24/im sample they are 3.5 and 1.7. This implies that the intermediate dust 
temperature component is more prominent in PAH dominated sources than in absorbed 
sources. For absorbed soures, if they contain a luminous AGN, the spectra can be flattened 
by having a more significant hot dust component to increase continuum emissivity at shorter 
wavelengths. 



Faint samples 



Finally, we compare our combined sample to those sources that are faint at observed frame 
24/im (f'24 < 2mJy). This is a difficult comparison to make as the IRS spectra of these 
sources are usually of low signal to noise, making detailed comparisons difficult. We there- 
fore make two adjustments. First, as most of our sample show PAH features, we restrict 
the co mparison to those 24/im sour ces that also show PAH features by using the compila- 



tion in 



Weedman and Houckl (120081 ). This compilation includes new spec tral measurements 



Branc 


et al.lbood: 


Weedman et al.l2006b: 


2006a 


Yan et al. 


2007; 


Pope et al. 


2008) 



Farrah et al.l2008; Houcket al 2005; Weedman et al 



Houck et al 



2007 
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used for faint sources at a variety of redshifts - the peak luminosity of the 7.7/xm PAH feature^ 
- and study how this peak luminosity evolves with redshift. The flux dens i ties f„ (7.7/im) and 
luminosities z/L I/ (7.7 y um) for our s ample are in Table El o r in brand et al. ( 2008 ). while those 
for the faint 24/im samples are in I Weedman and Houckl (120081 ) . 



We plot these 7.7/zm PA H luminosities against redsh ift in Figure El We also include 
the 24/xm br i ght so urces from lWeedman and Houckl (120081 ). and low redshift ULIRGs from 
Spoon et al.l (120071 ) . Two important results can be seen. First, the 70/im selected star- 
bursts are among the most luminous known from any infrared se lected samples. T hey are 
more luminous, on average, than both low redshift ULIRGs or the Yan et al. J2007h sources 
at comparable redshifts, and approach the luminosities of the 24/im and sub-mm selected 
sources at z > 1.5. This is expected - the high redshifts of our sample mean we are probing 
a greater volume and hence can find more luminous sources than local examples, and the 
additional demand of a 70/im detection means our sources will be more luminous, on average, 
than sources with just a 24/im detection at similar redshifts. 

Second, they reside in a redshift range, 0.5 < z < 1.5, where other selection methods 
turn up relatively few sources. The faint 24/im samples, which have comparable 24/im fluxes 
to our sample but are not detected at 70/im, span a significantly broader redshift range of 
0-5 ^-2^3, with the majority lying at z > 1.5. It seems therefore that a 70/im selection 
of ~20mJy, together with a faint optical and 24/im flux cut, serves to select sources almost 
entirely in the cruci al redshift range 0.5 < z < 1-5. This is the redshift range, for example, in 
which the results of lLe Floc'h et al.l (120051 ) show steeper evolution of the IR-luminous galaxy 
population than is shown in Figure [6j We conclude that a 70/im selection, in combination 
with selections at mid-IR and far-IR/sub-mm wavelengths, is vital to measure adequately 
the luminosity evolution of luminous starburst galaxies over < z < 2.5. 



5. Summary 



Sixteen spectra have been obtained with the Spitzer IRS of extragalactic sources in the 
SWIRE Lockman Hole field having f y (70/im) > 19 mJy, including 12 sources with optical 
magnitudes m r > 23. Results are combined with the sample of 11 s ources with i v {7Qu m) 
> 30 mJy from the NOAO Deep Wide-Field Survey region in Bootes (IBrand et al.l 120081 ) to 
consider the nature of the 70/tm population. 



As we are using the p eak luminosity o f this feature, ra t her th an its integrated flux, our luminosities 
differ from those quoted in Yan et al. ( 2007 ) and Pope et al. ( 20081) 
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The 70/iin sources are characterised either by strong PAH features or by strong silicate 
absorption features with weak or absent PAHs. Ten of the 16 objects show prominent PAHs; 
four show strong silicate absorption, and two have no discernable spectral features. The 
continuum luminosities (measured by z/L I/ (15/xm) in ergs s _1 ) span 43.8 < log z/Lv(15/zm) < 
46.3 with the absorbed sources having higher luminosities than the PAH dominated sources. 
Compared to sources that are bright at 24/xm (fj,(24/im) > lOmJy), the 70/zm sources have 
steeper rest frame mid-IR continua and higher luminosities. 

The 70/zm sources with strong PAH features are among the most luminous starbursts 
seen at any redshift. Furthermore, these sources effectively fill the redshift range 0.5 < z 
< 1.5 where previous selection methods using a 24yum flux of ~ lmJy but without a 70/iin 
detection have found few sources. This result demonstrates that selection of sources at 
70yum to fainter flux limits will provide crucial samples for determining the evolution of star 
formation with redshift. 

We thank the referee for a very helpful report. This work is based on observations 
made with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, 
California Institute of Technology under NASA contract 1407. Support for this work by 
the IRS GTO team at Cornell University was provided by NASA through Contract Number 
1257184 issued by JPL/Caltech. Support for the Spitzer Space Telescope Legacy Science 
Program, was provided by NASA through an award issued by the Jet Propulsion Laboratory, 
California Institute of Technology under NASA contract 1407. The research described in 
this paper was carried out, in part, by the Jet Propulsion Laboratory, California Institute of 
Technology, and was sponsored by the National Aeronautics and Space Administration. This 
research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated 
by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the 
National Aeronautics and Space Administration. DF thanks the Science and Technologies 
Facilities Council for support via an Advanced Fellowship. 
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Table 1. 


Observations Lo 


g 




ID 




Object 


Modules a 


Exposure times 
(s) 


AOR Key b 


1 


SWIRE4 


J103637.18+584217.0 


SL1/LL2/LL1 


240/480/480 


17410304 


2 


SWIRE4 


J103752. 14+575048.6 


SL1/LL2/LL1 


120/240/240 


17413888 


3 


SWIRE4 


J103946.28+582750.7 


SL1/LL2 


300/720 


17413120 


4 


SWIRE4 


J104057.84+565238.9 


SL1/LL2 


300/720 


17411840 


5 


SWIRE4 


J104117.93+595822.9 


SL1/LL2 


240/480 


17410560 


6 


SWIRE4 


J104439.45+582958.5 


SL1/LL2 


240/600 


17412608 


7 


SWIRE4 


J104547.09+594251.5 


SL1/LL2/LL1 


240/600/600 


17413376 


8 


SWIRE4 


J104827.68+575623.0 


SL1/LL2 


300/720 


17410048 


9 


SWIRE4 


J104830.58+591810.2 


SL1/LL2 


240/600 


17413652 


10 


SWIRE4 


J104847. 15+572337.6 


SL1/LL2/LL1 


180/360/360 


17412864 


11 


SWIRE4 


J105252.90+562135.4 


SL1/LL2/LL1 


240/600/600 


17410816 


12 


SWIRE4 


J105404.32+563845.6 


SL1/LL2/LL1 


120/240/240 


17412352 


13 


SWIRE4 


J105432.71+575245.6 


SL1/LL2 


300/720 


17409792 


14 


SWIRE4 


J105509.00+584934.3 


SL1/LL2 


300/480/720 


17412096 


15 


SWIRE4 


J105840.62+582124.7 


SL1/LL2/LL1 


240/600/600 


17414144 


16 


SWIRE4 


J105943.83+572524.9 


SL1/LL2/LL1 


240/600/600 


17411328 



a Modules of the Spitzer Infrared Spectrograph used to observe source. 

b Astronomical Observation Request number. Further details can be found by refer- 
encing these numbers within the Leopard software, available from the Spitzer Science 
Center. 



-19- 



Table 2. Basic Data 



ID z phot a z irs h m r IRAC Fluxes (>Jy) MIPS Fluxes (mJy) 















3. 6/i in 


4. 5/i in 


5. 8/i in 


8/i HI 


24/im 


70/im 


160/i hi 


1 






0.97 






21.4 


25.7 


48.8 


197.1 


2.17 


34.7 




2 


(0 


47) 


1.55: 


24 


06 


86.4 


248.2 


533.8 


1248.2 


3.98 


19.9 




3 


1 


30 


0.90 


23 


36 


28.0 


25.2 






1.05 


22.8 




4 






0.93 


22 


53 


79.0 


57.6 


47.1 




1.05 


24.2 




5 






0.65 


20 


83 


88.9 


77.8 


92.4 


123.5 


1.46 


32.9 


77.4 


6 






0.68 


21 


49 


46.1 


49.5 


59.0 


139.1 


1.21 


22.0 




7 





52 


0.39 


20 


21 


115.9 


105.5 


129.7 


319.3 


1.77 


20.2 




8 


1 


33 


0.86:: 






40.7 


41.0 


33.3 


138.0 


0.98 


37.4 




9 


1 


14 


0.94 


23 


28 


146.9 


122.8 


114.4 


114.4 


1.62 


20.6 


125.0 


10 


(1 


22) 


1.47: 


23 


74 


154.0 


286.1 


455.7 


799.1 


2.62 


23.2 


124.0 


11 






1.29:: 


20 


39 


236.6 


165.1 


164.9 


176.2 


1.53 


30.2 


71.3 


12 






1.72: 


23 


29 


23.4 


27.2 




102.0 


4.25 


26.4 




13 






1.02 


23 


29 


72.4 


66.6 


64.9 


113.0 


1.18 


37.0 


75.9 


14 






0.88 






68.4 


63.9 


67.9 


120.4 


0.97 


24.1 




15 


1 


06 


0.89 


23 


85 


112.3 


83.6 


89.4 


92.6 


1.50 


19.3 




16 


1 


18 


0.80 


23 


16 


170.8 


180.0 


218.2 


310.3 


1.96 


30.7 


119.8 



Note. — IRAC fluxes have errors of 5%; MIPS 24/im fluxes have errors of 10%, and 
MIPS 70/im and 160/iin fluxes have errors of 20%. 



a Photometric redshift, derived using the code of iRowan- Robinson et all (120081 ). Red- 
shifts in brackets are based on two optical bands and are not considered reliable. 

b Redshift derived from the IRS spectrum. These are accurate to Az = 0.02, except 
sources noted by ':', which are accurate to Az ~ 0.2, and those noted by'::' which have 
A z ~ 0.4 and should be regarded with caution. 



Table 3. Spectral Measurements 



ID 


PAH 6.2/im 


f7.7 a 


uL 7 . 7 b 


PAH 11.2/zm 


Ssil 


l24/fl5 C 


fl5 d 


log[^Li5] c 


SFR f 






Flux 


EW 


mjy 


ergs s — 1 


Flux 


EW 






mjy 


ergs s 


M yr _1 


logfLr^ ) 

o V W / 


1 


4.20±0.65 


0.18±0.06 


<3.1 


<45.41 


<2.0 


<0.3 


2.30±1.00 


(2.4) 


4.8 


45.31 


<690 


12.41 


2 


<1.10 


<0.15 


<5.2 


<46.10 


<4.0 


<0.2 


0.39±0.11 


(2.4) 


6.8 


45.92 


<3400 


13.01 


3 


1.55±0.41 


0.15±0.10 


0.75 


44.79 






<1.61 


(3.77) 


(1.76) 


44.87 


170 


12.23 


4 


4.31±0.80 


0.82±0.22 


1.05 


44.97 






<2.00 


(3.88) 


(1.86) 


44.92 


250 


12.28 


5 


4.29±0.83 


0.47±0.15 


1.55 


44.82 


<4.5 


<0.5 


0.87±0.80 


(3.85) 


(1.63) 


44.55 


180 


12.06 


6 


2.80±0.60 


0.31±0.20 


1.15 


44.73 


3.78±0.81 


0.68±0.35 


1.18±0.90 


(3.50) 


(1.41) 


44.53 


110 


12.02 


7 


1.79±0.25 


0.17±0.08 


<1.1 


<44.21 


4.10±1.00 


0.55±0.25 


1.25±0.70 


2.4 


1.4 


44.02 


<44 


11.48 


8 






0.7 


44.72 








(4.82) 


(1.68) 


44.81 


110 


12.30 


9 


9.04±1.18 


1.00±0.22 


3.4 


45.49 








(2.99) 


(2.61) 


45.08 


830 


12.55 


10 


<6.00 


<0.50 


<4.7 


<46.01 


<7.0 


<0.5 


0.81±0.25 


(5.2) 


3.3 


45.55 


<2700 


13.15 


11 
























12.80 


12 


<5.60 


<0.30 


<6.8 


<46.30 


<8.5 


<0.76 


1.32±0.24 


(2.2) 


(10.7) 


46.20 


<5400 


13.11 


13 


3.62±0.38 


0.58±0.09 


1.65 


45.24 








(5.22) 


(2.17) 


45.07 


470 


12.50 


14 


1.57±0.45 


0.35±0.14 


2.1 


45.22 






0.80±0.60 


(3.97) 


(1.63) 


44.82 


450 


12.35 


15 


7.19±0.31 


1.27±0.18 


2.8 


45.36 


<9.00 


<2.3 


0.69±0.35 


(4.57) 


1.50 


44.79 


620 


12.36 


16 


5.20±0.81 


0.36±0.10 


2.3 


45.18 


4.79±1.66 


1.24±1.00 


1.94±1.00 


(3.70) 


2.40 


44.90 


410 


12.48 



O 



Note. — For the 6.2/xm and 11.2/xm PAH features, fluxes are in units of 10 — 22 W cm -2 and rest frame equivalent widths in /jm. 
a Observed frame flux density at the peak of the 7.7/^m PAH feature. Error is approximately 10%. 
b Rest frame PAH luminosity determined from the peak f„(7.7/xm). 

c Rest-frame f^(24/im)/f l ,(15/im) continuum slope. Measurements are made from the IRS spectra if both rest-frame wavelengths are seen, and 
have a ~ 10% error; otherwise, fluxes at one or both wavelengths are estimated via interpolation from a power law with a slope determined from 
comparison of observed f„(24/im) and f„(70/xm). Ratios determined using such interpolations are in parentheses, and have errors of ~ 20%. 

d Rest-frame 15^tm flux density. Sources for which this is interpolated via a power law are in parentheses. 

c Rest frame 15^tm continuum luminosity vL v (l5fim). Sources for which this luminosity is interpolated via a power law are in parentheses. 
f Star formation rate, determined from Equation [T] 

g Rest-frame 1-1000/^m luminosity derived from the SED fits described in £|3.3I The error on the luminosities is approximately 25% in all cases, 
and does not include errors arising from uncertainties in redshift. 
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Fig. 1. — Individual spectra of sources with wavelength in /im and flux density in mJy, 
plotted in the rest-frame using the IRS redshifts in Table [2j The numbers in the top left of 
each panel are the ID number in Table Q] followed by the redshift. The vertical dotted lines 
mark the 6.2/zm, 7.7/zm, 8.6/xm, 11.2/im and 12.7/im PAH features. 




Fig. 2. — Individual spectra, continued from Figure [T] 
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Fig. 3. — Left panel: Redshift histogram, using the IRS redshifts. The soli d line is the 
distribution for our sample. The dotted line is the distribution for the sample of iBrand et al. 



(120081)- Right panel: Comparison between the IRS redshifts and the photometric redshifts 
from Rowan-Robinson et al. fboosh (Table [2]). The dotted lin es denote a deviation of 0.06 i n 
log(l+z), the boundary defined as a 'catastrophic failure' by iRowan- Robinson et al.l (120081 ). 
Open symbols (objects 2, 8 and 10) have an unreliable IRS and/or photometric redshift. 
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Fig. 4. — The 'Fork' diagram of ISpoon et al.l (120071 ). Tria ngles: Our 70 am so urces. Cir- 
cles: Local ULIRGs. Stars, squares: Classical starbursts (IBrandl et al.l 120061 ) and AGN 
(jWeedman et al.ll2005l ) with IR luminosities in the range 10 10 - lO n L , respectively. 
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Fig. 5. — Distribution of continuum luminosity (uL u (15fim) vs. rest rame continuum slope 
(f I/ (24/im)/f I/ (15/xm)). Error bars have been omitted for clarity, but are in most cases compa- 
rable to the symbo l sizes . Filled symbols are sources from our combined 70/zm sample (Table 



3 and Brand et al. 



(200 g)), while open symbols are sources from the flux limited, f iy (24/zm) > 



10 mJy sample of 
Triangles: PAH dominated systems. 



Weedman and Houckl (120091 ). Squares: sources without detectable PAHs. 
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Fig. 6. — Distribution of PAH luminosity uL u (7.7fim) with redshift. Error bars have been 
omitted for clarity, but are in most cases comparable to the symbol sizes. R ed Triangles 



PAH dominated sources from the combined 70 /im sample from Table 3 and iBrand et al. 



(hoosh. Black Circl es: Starburst component of l ow redshift ULIRGs flSpoon et al. 



Imanishi et al.ll2007l ). Blue Squares: Starbursts in lWeedman and Houckl (120081). Solid line: 



2007 



l umin osity evolution to z — 2.5 of the form (1 + z) 2 ' 5 , determined in IWeedman and Houck 
(120081 ) . Results show that 70/im PAH sources are among the most luminous starbursts in 
the Universe, and fill a crucial redshift regime. 



